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1 Introduction 



Charm mixing in the Standard Model is described by the parameters x = 2(M 2 — 
Mi)/(T 2 + r\) and y = (T 2 — ri)/(r 2 + Ti), where M 1>2 and r 12 are the masses 
<D and widths, respectively, of the neutral D meson CP eigenstates, D\ (CP-odd) and 

D 2 (CP-even): \Di }2 ) = (|P°) ± \D°))/\^2. Recently, several experiments have 
directly probed x and y [H EJ |3], IU [5j EJ [7], as well as the "rotated" parameters 
y' = y cos 5 — x sin S and x' = y sin 5 + x cos S [H [91 [10] , where —5 is the relative phase 
y—i between the doubly Cabibbo-suppressed D° — > i^ + 7r _ amplitude and the correspond- 

ing Cabibbo-favored P° K + n- amplitude: (K + n-\D°) / (K + n-\D°) = re' 15 . We 
O adopt a convention in which 5 corresponds to a strong phase, which vanishes in the 

SU(3) limit [TTj. In this article, we update an analysis [12] that gave the first direct 
determination of cos 5 using correlated production of D° and D° mesons in e + e" col- 
lisions produced at the Cornell Electron Storage Ring and collected with the CLEO-c 
detector. Here, we also present a first measurement of sin 5. 

At the ^(3770) resonance, the D°D° system is produced in an overall charge 
conjugation eigenstate with C = — 1. As a result of this quantum coherence, the 
exclusive D°D° decay rate to a final state {i, j}, where i and j label the final states 
of the two D mesons, is given by the square of the antisymmetric amplitude 



i,j)<xMl= (^|P )(j|P U )-(z|P U )(j|P ) +0(x 2 ,y 



,2> 



where the 0(x 2 , y 2 ) term represents a mixed amplitude. These exclusive rates depend 
on the mixing parameters and amplitude ratios defined above: 

T(i, j) = T(z, j) oc 1 + r 2 r 2 — 2ri cos 5i Vj cos 5j — 2r, L sin 5i Vj sin 5j (2) 
T(i, j) = r(z, j) oc r 2 + r 2 — 2r { cos <5j r 3 - cos 5j + 2r^ sin 5i Tj sin 5j. (3) 
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Inclusive rates are obtained by summing over exclusive rates and are trivially related 
to the branching fractions (£>) in uncorrelated decay: 

F(i,X) = F(hj) + =B i + Bi<xl + 2yr i cos6 i + rf . (4) 

j 

By comparing correlated and uncorrelated decay rates, we extract simultaneously the 
charm mixing and strong phase parameters, as well as the number of D°D pairs 
produced (TV) and the branching fractions of the reconstructed D° final states. 

Compared to our previous analysis [12], this update uses a dataset three times 
larger (818 pb~ x ), and we reconstruct more CP eigenstates, as well as semimuonic 
D° decays and modes that provide sensitivity to sin S and r. 

2 Experimental Technique 

We reconstruct the D meson final states listed in Table [TJ Following Ref. [13] . the 
Kg7T + 7r~ Dalitz plot is divided into eight bins according to the strong phase of the 
decay amplitude. Because CP eigenstates and semileptonic final states have known 
values of T\ and Si, they act as reference points for determining y and the strong phase 
S in the Ktt final state. Inclusive rates are derived from yields of single tags (ST), or 
individually reconstructed D° or D° candidates. For exclusive rates, we reconstruct 
D°D pairs or double tags (DT). 

Fully reconstructed modes are identified using two kinematic variables: the beam- 
constrained candidate mass M = yEfi /c 4 — p|, / c 2 and the energy difference AE = 
Ed — Eo, where Pd and Ed are the total momentum and energy of the D candidate, 
respectively, and E is the beam energy. We measure ST yields for all fully recon- 
structed modes except i^^7r + 7r _ . Decays with K° L mesons and neutrinos, which we 
do not detect directly, are identified with a partial reconstruction technique, where 
the presence of the undetected particle is inferred via conservation of energy and 
momentum. We form DT candidates from all combinations of modes in Table [TJ 



Type 



Reconstruction 



Final States 



S; 



Mixed-CP 
CP+ 
CP+ 
CP- 
CP- 
Semileptonic 



Full 
Full 

Partial 

Full 
Partial 
Partial 



, 7T + 7T", K° s n°TT° 



K-7T+, K+71 

K + K 

K° L 7T , K° L 7], K° L u 

K° s 7i°, K° sV , K° s u 
K° L ir°ir° 
K~{e + ,n + }v e ^, K + {e~,fi-}v 



From fit 
1 
1 
1 
1 




7T 
71 






Table 1: D final states reconstructed in this analysis. 
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both Cabibbo-favored and Cabibbo-suppressed, with at most one missing particle. 
In addition, we reconstruct {Kev e , K^tt }, which has two missing particles, using the 
technique described in Refs. [T4l [To]. 



3 Results 

We combine 261 yield measurements, along with estimates of efficiencies and back- 
ground contributions, in a least squares fit that determines 51 free parameters, in- 
cluding the five charm mixing and D — > Ktt amplitude parameters listed in Table |2j 
The other parameters, for which the fit results are not shown, are A/", amplitude ra- 
tios and phases for the 8 phase bins in Kg7r + ir~, and 21 branching fractions. The 
fit includes both statistical and systematic uncertainties on the input measurements. 
We separate the statistical and systematic uncertainties in Table [2] by repeating the 
fit with only statistical uncertainties on the inputs. We perform one fit with no ex- 
ternal inputs (Standard Fit) and another including the external measurements of y, 
x, r 2 , y', and x' 2 compiled in Refs. [TOJ [T7J (Extended Fit). In the Standard Fit, 
there is a sign ambiguity for sin S, which is resolved in the Extended Fit. Because 
we treat uncertainties on external measurements as systematic uncertainties, when a 
fit parameter is directly constrained by an external measurement, we quote only one 
uncertainty in Table [2j in these cases, the statistical uncertainty effectively vanishes. 

In the Standard Fit, the statistical uncertainties on y and cos 5 are roughly three 
times smaller than in our previous analysis [12J . The Extended Fit demonstrates that 
our measurements of cos 5 and sin 5 can be used to improve the uncertainty on y (by 
combining y and y' from other experiments) by approximately 10%, compared to the 
current world average found in Ref. [17]: y = 0.79 ± 0.13. 

Asymmetric uncertainties are determined from the posterior probability distribu- 
tion functions (PDFs) shown in Fig. [TJ These curves are obtained by re-minimizing 



Parameter Standard Fit Extended Fit 

y (%) 4.2 ±2.0 ± 1.0 0.636 ±0.114 

r 2 (%) 0.533 ±0.107 ±0.045 0.333 ± 0.008 

r q oi +0.22+0.07 +0 .19+0.00 

LOb0 u - oi -0.18-0.05 i - id -0.17-0.08 

sin S -0.01 ± 0.41 ± 0.04 0.56±g;|f±^;lJ 

x 2 (%) 0.06 ±0.23 ±0.11 0.0022 ± 0.0023 



Table 2: Results from the Standard Fit and the Extended Fit for all parameters 
except branching fractions. Uncertainties are statistical and systematic, respectively. 
In the Extended Fit, we quote only one uncertainty for y, r 2 , and x 2 , which are 
directly constrained by an external measurement. 
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the x 2 & t each point and computing C = e ^ x miJ/ 2 . We construct the PDFs for 5 
by scanning cos 5 and sin 5 under the constraint cos 2 5 + sin 2 5 = 1, which result in 
implied values for 5 of \5\ = (10i|lJ 3 )° for the Standard Fit and 5 = (18i#)° for the 
Extended Fit. 



3950913-009 3951913-010 




sin8 cosS sin5 cos5 



Figure 1: Likelihoods for the Standard Fit (left) and Extended Fit (right) including 
both statistical and systematic uncertainties for cos 5, sin 5, 5, and y. The two- 
dimensional likelihoods are shown as solid contours in increments of la, where a = 
\J Ax 2 - The dashed contour marks the physical boundary, where cos 2 5 + sin 2 5 = 1. 

4 Summary 

We present an updated analysis of charm mixing and the D — > Kir strong phase 
using quantum correlations in D°D decays at the ^(3770) resonance. These results 
are based on the full CLEO-c dataset, and they make a significant contribution to 
the world averages of mixing parameters. 
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